A Time-Like Naked Singularity 
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We construct a class of spherically symmetric collapse models in which a naked singularity may 
develop as the end state of collapse. The matter distribution considered has negative radial and 
tangential pressures, but the weak energy condition is obeyed throughout. The singularity forms 
at the center of the collapsing cloud and continues to be visible for a finite time. The duration of 
visibility depends on the nature of energy distribution. Hence the causal structure of the resulting 
singularity depends on the nature of the mass function chosen for the cloud. We present a general 
model in which the naked singularity formed is timelike, neither pointlike nor null. Our work 
represents a step toward clarifying the necessary conditions for the validity of the Cosmic Censorship 
Conjecture. 



PACS numbers: 04.20.Dw, 04.70.-s, 04.70.Bw 

The cosmic censorship conjecture (CCC) has been 
widely recognized as one of the most important open 
problems in gravitational physics today. This is because 
several important areas in the theory and applications 
of black hole physics crucially depend on CCC. Nev- 
ertheless, the CCC remains unproved and there exists 
no mathematically precise and definite statement for the 
CCC which one could try to prove (see e.g. [1-6] for some 
recent reviews, and references therein). 

For this reason a detailed study of dynamically devel- 
oping gravitational collapse models within the framework 
of general relativity becomes rather essential. The hope 
is that such a study may allow us to formulate a provable 
statement of the CCC, if it is correct in some form. Such 
investigations also help us to discard certain statements 
of the CCC which might sound plausible but for which 
there exist counter-examples which show that the CCC 
cannot be valid in such a form. They may even illustrate 
the physical conditions that give rise to naked singulari- 
ties (NS) or black holes (BH) as end states of a realistic 
gravitational collapse. So far, such dynamical collapse 
studies have focused largely on collapse models that cre- 
ate either BH or NS, depending on the nature of the ini- 
tial profiles of density, pressure, and velocity from which 
the collapse develops. In many of these cases, when a 
NS develops, it is located at the center of the spherically 
symmetric cloud (a central singularity, see e.g. [5-7]). In 
that case, there will exist families of non-spacelike future 
directed geodesies, which will be accessible to distant ob- 
servers in the future, and which will terminate at the 
singularity in the past, thus making it visible in princi- 
ple. This is opposed to the BH case where the apparent 
horizon forms early enough to cover all of the singular- 



ity, with no portion of it remaining visible to outside 
observers. 

If we require the pressure to be positive then the " cen- 
tral" singularity, if it is naked, corresponds to a singular- 
ity along a visible null line. The remainder of the singu- 
larity is spacelike and covered by a horizon. In this paper, 
however, we permit the pressure to be negative and ex- 
amine the structure of the singularity. We construct an 
explicit solution in which the singularity may be time- 
like. It may even change its character, being timclike 
along a certain region and, after being visible for a finite 
time, turning spacelike and being covered. The collapsing 
matter is described by a particularly chosen matter field 
that satisfies the weak energy condition although the ra- 
dial and tangential pressures are negative and unequal. 
While what is presented here is a specific construction 
of a class of collapse models, involving somewhat special 
choices, we make sure that physical reasonability condi- 
tions such as the energy conditions and the regularity of 
the initial data at the initial surface are respected. 

The spherically symmetric metric in a general form can 
be written as 

ds 2 = -e 2 ^dt 2 + e 2 ^dr 2 +R 2 (t,r)dn 2 , (1) 

where dfl 2 is the line element on a two-sphere. Choosing 
the comoving frame, the stress-energy tensor for a general 
(type I) matter field is given in a diagonal form as 

Tl = -p; r; = Pr ; T e B = T* = p B , (2) 

where p, p r and pg are the energy density, and the radial 
and tangential pressures respectively. We assume that 
the matter field satisfies the weak energy condition, that 



2 



is, the energy density as measured by any local observer 
is non- negative and, for any timelike vector V, 



TikVW* > 



which amounts to, 



P > 0; p + Pr > 0; p 



> 0. 



(3) 



(4) 



The initial data consists of three metric functions, the 
energy density, and the radial and tangential pressures 
at the initial time t = t,. This is given in terms of 
six arbitrary functions of the radial coordinate, viz., 
v(U,r) = vo(r), ip{U,r) = ip (r), R(U,r) = r, p(U,r) = 
Po(r), Pr (*»>!") = Pro(r), Pe(U,r) = Peo(r), where, us- 
ing the scaling freedom for the radial co-ordinate r we 
have chosen R(ti,r) = r at the initial epoch. The dy- 
namic evolution of the initial data is then determined by 
the Einstein equations, which for the metric Q become 
(8ttG = c = 1), 



F' 
R 2 R' 



2{pg - Pr ) R' 
p+Pr R 



R 2 R 



Pr 



-2R' + R'% + R^=0 



G - H=l -R> 



(5) 



(6) 



(7) 



(8) 



where F = F(t, r) is an arbitrary function. In spherically 
symmetric spacetimes F(t,r) is interpreted as the mass 
function, with F > 0. In order to preserve the regular- 
ity of the initial data we must also require Ffa, 0) = 0, 
i.e., the mass function should vanish at the center of the 
cloud. The functions G and H are defined as G(t, r) = 



e- 2 v(R'f and H(t,r) 



All the initial data above are not mutually indepen- 
dent: from equation © wc find that the function Vo{r) 
is determined in terms of rest of the initial data. Also, 
by rescaling of the radial coordinate r, the number of 
independent initial data functions reduces to four. We 
then have a total of five field equations with seven un- 
knowns, p, p r , pe, ip, v, R, and F, giving us the freedom 
of choice of two free functions. Selection of these func- 
tions, subject to the given initial data and weak energy 
condition, determines the matter distribution and metric 
of the space-time and thus leads to a particular collapse 
evolution of the initial data. At this point it is convenient 
to introduce a scaling variable v(t, r), defined as 



R(t,r)=rv(t,r), (9) 



where, 



the last condition being necessary for a collapse. Let us 
consider the following choice of the allowed free functions, 
F(t, r) and u(t, r), 



F(t,r) = r 3 M{r)v, 



(11) 



where r 3 A4(r) is a suitably differentiable and monotoni- 
cally non-decreasing function, and 



v(t,r) = vq(R). 



(12) 



The function M may be expanded in a Taylor series 
about r = 0, 



M(r) = M Q + M 2 r 2 + M 3 r 3 + 
Then, from equation JSJ, we have 

3Mv + r[M, r v + Mv'] 



P = 



and 



Pr 



v 2 \v + rv'] 



M(r) 



(13) 



(14) 



(15) 



The above choice of mass function therefore imples that 
the radial pressure is negative. The weak energy condi- 
tion, however, does hold. If R' = v + rv' and F' are 
both positive in the equation l|14|l. then we clearly p > 0. 
Again, for p + p r > at all epochs, it must be true that 

(2M Q + AM 2 r 2 + 5M 3 r 3 + ••>> 0. (16) 

But, because v > 0, it follows that if the condition 
p + p r > is satisfied at the initial epoch, it is satisfied 
throughout the evolution. Finally, from equation JG|), 



1 



2 M 



P + Pe = ^(P + Pr) [1 + Rv (R),r] + > (17) 

if [1 + Rvo(R) t ii] > for all epochs. This provides a 
necessary condition for the weak energy condition to be 
satisfied. 

At the initial epoch we then have 



and 



Po (r) = 3M + 5M 2 r 2 + 6M 3 r 3 + 



Pr {r) =-[M + M 2 r 2 + M 3 r 3 + 



(18) 



(19) 



v(U,r) = l, v(t s (r) i r) = and v < 0, (10) 



At the initial epoch, the radial and the tangential pres- 
sures must be equal at the center and all the pressure 
gradients must vanish. It follows that the initial tangen- 
tial pressure must have the form 

Pe (r) = -[M Q + P e 2 r 2 + P e 3 r 3 + ■■■}■ (20) 

Hence from equation © we see that Vo(r) becomes, 

Mr) = a 2 r 2 + a 3 r 3 + ■■■ (21) 
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where 



a-2 



Mi 



a 3 



Pe 3 + M 3 
2M 



2M 

and, from equation (|12fl . 

v(t, r) = v (R) = a 2 R 2 + a 3 R 3 



(22) 



The dynamic evolution of pe(t, r) is obtained by inserting 
equation l|12|) in equation© and simplifying to get, 



I \ , R Pr 

Pe(r,v) = p r + — 



1 



v (R) tR R(p + p r ) 



(23) 



There therefore exists an e ball around the central shell 
for which pg — p r and the perfect fluid equation of state 
is valid. 

Using equation 11211 in equation Q, we get, 



G(t,r) = b(r)e 2 " a(R \ 



(24) 



where b(r) is another arbitrary function of r. In corre- 
sponding dust models, we can write b(r) = 1 + r 2 bo(r), 
where bo(r) is the energy distribution function of the col- 
lapsing shells. Thus, the metric Q becomes, 



dg 2 = _ e 2 ( a 2 R- + -) dt 2 



R >2 e -2(g 2 R 2 +-) dr 2 

l + r 2 b n (r) 



(25) 

and is valid for small values of r, for all epochs, i.e. for 
all values of v(r,t), till the singularity. 

Solving the equation of motion (JSJ) we find that 

v = -e 2uo{rv) ^v 2 {2a 2 + 2a 3 rv ■■■)+ b (r)e 2 ^irv) + M(r), 

which may be integrated to obtain 



t(v,r) 



e -2» (rv) dv 



Iv y / v 2 (2a 2 + 2a 3 rv ■■■)+ b Q {r)e 2v ^ rv ) + M(r) 

(27) 

We note that the radial coordinate r is treated as a con- 
stant in the above equation, which gives the time taken 
for a shell labeled r, to reach a later epoch v in collapse 
from the initial epoch v = 1. 

It is clear that an explicit solution of the above integral 
will give a closed form solution of the form t = f(v,r) or, 
inversely, v = g(t, r), which will then determine the met- 
ric function R = rg(t, r) thereby giving an exact solution 
for the metric l|25|l . Unfortunately, the integral cannot be 
expressed in closed form and we make a Taylor expansion 
of the integral about the center of the cloud. 



t(v, r) = t(y, 0) + rX(v) + 0{r 2 ) 
where the function X{y) is given by, 



X{v) 



dv- 



2v 4 a 3 + h 



v /6 (0) + 2w 2 a 2 +7W(0) 



(28) 



(29) 



If a closed form solution of R exists up to the first ap- 
proximation, it will be of the form R = rX^ 1 [{t{v, r) — 
t(v,0))/r]. Therefore, by expanding as above we are ac- 
tually solving for R and so for the metric (|25|l to the 
first approximation, although we do not write it in closed 
form. This is because it is only the sign of ^(0) that de- 
termines the final end state of the collapse, which is the 
issue of interest here. 

The time taken for the central shell at r = to reach 
the singularity, t s (0), is given by 



t a (0) = 



dv 



o ^v 2 (2a2) + b (0)+M(0) 



(30) 



The time taken for the other shells (r ^ 0) to reach the 
singularity, t s (r), can be given as, 



t s (r) =t s (0) + rX{0)+O(r 2 ), 
where the function X(0) is given by, 



X(0) 



1 



dv- 



2v 4 a 3 + bi 



v /6 (0) + 2v 2 a 2 +M(0) 



(31) 



(32) 



We see that by suitable choice of the coefficients of initial 
density, pressure and energy profiles we can make ^(0) 
positive or negative. Furthermore, at the singularity, for 
a constant v surface we have 



lim v 



VMr) +M(r)X(Q) + 0(r 



(33) 



and, because we have expressions for v and v, near the 
central singularity, we can in principle calculate v(r, t) in 
the neighborhood of the central singularity. This solves 
the system of Einstein equations. 

The apparent horizon is the boundary of the trapped 
region of the space-time and is given by R/F = 1. If 
the neighborhood of the center gets trapped earlier than 
the singularity, then it will be covered and a black hole 
will be the final state of the collapse. Otherwise, the sin- 
gularity can be naked with non-spacelike future directed 
trajectories escaping from it to outside observers. Using 
(|11|1 . we find that the apparent horizon is just the surface 



r 2 M(r) = 1. 



(34) 



Therefore, if r 2 M.(rb) < 1 there will be no trapped sur- 
faces in spacetime, where rt, is the radial coordinate of 
the boundary of the cloud. 

It is simplest to examine the nature of the singularity 
by noting that it occurs at R = 0. This implies that at 
the singularity, 

R' 2 

ds 2 = [exp(2?A) - exp(2^)-^-]dr 2 . (35) 

If the right hand side is negative the singularity is time- 
like. Therefore, for a timelike singularity, G — H > 0, or 



1 -r 2 M{r) > 0. 



(36) 
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But, because the function r 2 A4(r) is monotonically non- 
decreasing, it follows that the singularity is timelike near 
r = 0, becomes null at r 2 M(r) — 1 and finally spacelike 
when r 2 M(r) > 1. 

This result may also be obtained explicitly near the 
center by examining the outgoing null geodesies. To see 
this specifically at r = 0, the outgoing radial null lines 
are given by 



then, for all values of r for which 



dt 
dr 



which, at the singularity, corresponds to 
= lim 



rv 



all Vl + r^r) 



(37) 



(38) 



In order to find the existence or otherwise of an outgoing 
null geodesic from the singularity we substitute the value 
of v' at the singularity in the above equation to obtain 



null 



If 



ry/bg(r) + M(r) 
VI + r 2 b {r) 



(39) 



(40) 



yfro(r) + M(r) 
^l + r 2 b {r) 



< 1 



(41) 



or 1 — r 2 M(r) > , there will be a visible outgoing null 
geodesic leaving the singularity. 

The model discussed here is based on a special choice 
of the mass function. The fluids described by this choice 
are not intended to describe an actual physical system: 
both the radial and the tangential pressures are negative. 
Nevertheless the initial data satisfy appropriate regular- 
ity conditions and the weak energy condition is main- 
tained throughout. Therefore, the system can serve as a 
guide to further clarifying and precisely formulating the 
CCC. 

Understanding what is possible in dynamically devel- 
oping collapse models is necessary to arrive at a plausi- 
ble concrete statement of the CCC. We have an example 
that shows that even in spherically symmetric collapse, 
the naked singularity need not always be either point- 
like or null, but can have an interesting causal structure, 
including being timclikc. 
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